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manager. 
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INVESTIGATION  OF  THE  INTERRELATIONSHIP  BETWEEN 
DIRECT,  DIFFUSE,  AND  TOTAL  SOLAR  RADIATION 


1 INTRODUCTION 


Background 

Solar  radiation  is  commonly  measured  by  its  total  intensity  on  a 
horizontal  surface.  Studies  of  solar  heat  gain  in  buildings  and  the 
performance  of  solar  collectors  require  information  about  total  radi- 
ation intensities  on  tilted  surfaces.  Horizontal  radiation  is  trans- 
formed to  tilted  surface  radiation  by  decomposing  the  total  horizontal 
radiation  into  its  direct  and  diffuse  components  and  then  recombining 
them  in  proper  proportions.  Direct  radiation  intensity  on  a surface  is 
given  by  the  direct  intensity  normal  to  the  sun's  rays  times  the  cosine 
of  the  angle  of  incidence  of  the  sun's  rays  on  the  surface.  Diffuse 
radiation  intensity  on  a surface  is  given  by  horizontal  diffuse  in- 
tensity times  the  radiation  view  factor  between  the  surface  and  the  sky. 

When  this  study  was  begun,  the  only  available  algorithm  for  com- 
puting diffuse  and  beam  radiation  from  the  total  horizontal  radiation 
was  based  on  daily  totals  of  the  three  components  listed  above.  Several 
computer  simulations  were  in  use  or  being  developed  which  computed 
hourly  building  loads  or  solar  performance.  These  programs  required 
hourly  values  of  beam  and  diffuse  radiation,  and  only  hourly  horizontal 
radiation  values  were  available.  It  was  not  known  how  accurately  the 
daily  algorithm  was  for  computing  hourly  diffuse  and  beam  radiation. 


Objective 

The  objective  of  this  research  was  to  determine  a simple  cor- 
relation between  total  radiation  on  a horizontal  surface  and  its  direct 
and  diffuse  components. 


Approach 

Using  two  solar  pyranometers,  researchers  measured  total  and  dif- 
fuse radiation  on  a horizontal  surface.  Data  collected  during  different 
seasons  showed  the  effects  of  daily  and  seasonal  changes  in  the  sun's 
position  and  of  seasonal  weather  changes.  These  data  were  recorded  by 
an  automatic  digital  system  which  created  data  tapes  suitable  for  com- 
puter analysis.  Several  regression  algorithms  were  used  to  find  which 
one  would  best  describe  the  measured  data. 
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Organization  of  Report 

Chapter  2 describes  the  equipment  and  the  procedure  used  to  measure 
the  horizontal  and  diffuse  radiation.  Chapter  3 presents  the  results  in 
terms  of  the  data  points  measured  and  the  calculated  correlation.  Chap- 
ter 4 analyzes  the  results  by  comparing  the  correlation  to  published 
correlations.  Chapter  5 gives  conclusions  and  recomnendations. 


2 EXPERIMENTAL  PROCEDURE 


Equipment  and  Its  Use 

Instrumentation  for  this  project  was  incorporated  into  instrumen- 
tation that  was  already  being  used  in  the  existing  solar  test  facility. 
The  solar  facility  had  one  pyranometer  for  measuring  total  radiation  on 
a horizontal  surface,  as  well  as  a system  for  automatically  recording 
signals  from  sensors  placed  throughout  the  building.  A second  pyr- 
anometer and  a shadow  band  were  purchased  to  measure  the  diffuse  radi- 
ation on  a horizontal  surface  by  the  shadow  band  blocking  the  beam  radi- 
ation. The  band  was  oriented  so  that  it  was  between  the  sun  and  the 
pyranometer  as  the  sun  moved  across  the  sky  during  the  day.  The  posi- 
tion of  the  shadow  band  was  adjusted  every  few  days  by  moving  the  shadow 
band  until  the  pyranometer  was  centered  in  its  shadow.  These  adjust- 
ments were  made  to  account  for  the  changes  in  solar  position  between 
summer  and  winter. 

A mounting  platform  was  installed  at  the  peak  of  the  solar  facility 
next  to  the  original  pyranometer.  This  location  provided  a clear  view 
to  the  horizon  in  all  directions.  The  new  pyranometer  was  connected  to 
the  solar  facility  data  recording  system. 

The  pyranometers  were  calibrated;  the  new  pyranometer  was  the  stan- 
dard, since  its  calibration  was  at  least  2 years  more  recent  than  that 
of  the  old  one.  The  voltage  signals  from  each  pyranometer  are  fed 
through  amplifiers  to  the  solar  facility's  data-recording  equipment. 

This  equipment  integrates  and  records  these  signals  and  38  other  signals 
every  hour.  Each  data  cassette  holds  approximately  1 week's  data. 
Because  of  limited  access  to  the  transfer  facilities,  several  cassettes 
were  accumulated  before  their  data  were  transferred  to  a large  tape. 

This  tape  was  then  transferred  to  a computer  facility. 


Procedure 


A program  was  written  to  compute  the  solar  position  and  the  value 
of  the  solar  constant  (intensity  of  the  sun's  direct  radiation  outside 
the  earth’s  atmosphere)  for  any  hour  of  any  day  of  the  year.  These 
values  were  combined  with  the  measured  horizontal  and  diffuse  radiation 
values  onto  a single  data  file.  This  file  contained  all  the  data  neces- 
sary to  analyze  and  evaluate  an  estimating  algorithm.  Another  program 
was  written  that  would  use  plotting  routines  to  generate  simple  plots  of 
the  data  points.  Other  programs  were  written  that  would  do  regressions 
against  different  algorithm  types  to  determine  the  best  one  to  use. 
Several  simple  polynomial  forms  were  tried,  because  polynomials  are  easy 
to  implement  in  simulation  programs. 
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The  data-gatheri no  technique  twice  allowed  errors  to  go  unnoticed; 
this  resulted  in  the  loss  of  several  months'  data.  First,  an  inter- 
mitter.-  amplifier  failure  was  not  observed  until  September  1977,  and 
there  was  no  spare  amplifier  available.  Return  to  the  manufacturer, 
repair,  and  reinstallation  took  until  the  end  of  January  1978.  Second, 
an  error  in  the  setting  of  the  cassette  recording  speed  during  April  de- 
stroyed the  next  3 months'  data. 

Experimental  results  are  presented  in  terms  of  the  following  quan- 
tities: 

Z = the  angle  between  the  sun  and  the  zenith 

H = the  total  radiation  measured  on  a horizontal  surface 

B = the  direct  (beam)  radiation  measured  normal  to  the 
direction  of  the  sun's  rays 

D = the  diffuse  radiation  measured  on  a horizontal  surface 

BO  = the  theoretical  beam  radiation  on  a surface  outside  the 

earth's  atmosphere 

HO  = the  theoretical  radiation  on  a horizontal  curface  outside 
the  earth's  atmosphere 

These  values  are  related  by  the  expressions 


H = B * cos(Z)  + D 

[Eq  1] 

HO  = BO  * cos(Z) . 

[Eq  2] 

1; 


Figures  1 through  4 show  the  experimental  data,  giving  the  ratio 
B/BC  plotted  as  a function  of  H/HO.  Figure  1 is  for  the  period  from 
February  1 through  March  5,  while  Figure  2 is  for  March  6 through  April 
6.  Figure  3 is  for  August  3 through  August  25,  and  Figure  4 combines 
the  data  from  Figures  1 through  3. 

Attempts  to  fit  the  data  with  a simple  polynomial  (of  the  form  Y 
= A°  + AiX  + A2X2  + A3X3  +...)  did  not  produce  satisfactory  results.  A 
simple  fit  was  made  by  computing  the  average  value  of  B/BO  for  each  10 
percent  increment  of  H/HO.  Table  1 provides  the  values  found  by  this 
fit.  Curves  created  by  interpolation  between  these  values  are  shown  in 
Figures  1 through  4. 
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L\  ANALYSIS 


Compari sons 

Several  current  computer  simulations  use  a correlation  which  has 
been  presented  in  equation  form  in  the  TRNSYS  solar  simulation 
program.  1 The  correlation,  which  is  based  on  work  by  Liu  and  Jordan,2 
is  predicted  on  daily  values  of  horizontal  and  diffuse  radiation  and  is 
expressed  as  a third-order  polynomial.  Several  other  radiation  cor- 
relations have  been  developed  recently.  One  was  used  in  a recent  com- 
parison of  several  solar  system  simulation  programs.3  Historical  hourly 
data  from  many  sites  were  used.  The  relationship  between  direct  and 
horizontal  radiation  is  expressed  by  a simple  linear  function.  A third 
method  was  developed  for  the  National  Climatic  Center  as  part  of  an 
effort  to  produce  standard  weather  data  tapes  called  SOLMET  tapes.4 
This  method  collected  data  from  more  than  20,000  points,  using  carefully 
calibrated  instruments  at  four  locations  in  the  United  States  having 
very  different  climates.  The  SOLMET  correlation  uses  the  averaging  pro- 
cess for  every  10  percent  increment  of  H/HO. 

Figure  5 shows  three  correlations  and  the  correlation  obtained 
during  this  study.  This  study  is  in  best  agreement  with  the  SOLMET  cor- 
relation, but  it  generally  predicts  less  beam  radiation.  Figure  6 com- 
pares the  correlations  for  the  time  periods  shown  in  Figures  1 through  3 
to  SOLMET.  This  illustration  shows  good  agreement  among  the  March, 
August,  and  SOLMET  data.  The  February  curve  is  significantly  lower  and 
shows  some  values  of  H/HO  which  are  much  higher  than  the  normal  range. 
Figures  1 through  3 show  a small  problem  with  the  averaging  algorithm  in 
the  .80  to  .90  H/HO  increment.  The  data  points  are  not  evenly  spread 
across  the  entire  range,  occurring  mostly  at  the  low  end.  This  causes 
the  value  of  B/BO  to  be  too  low  in  this  range.  The  SOLMET  correlation 
shows  an  even  lower  value  in  this  increment,  because  it  had  relatively 
fewer  data  points  above  the  .82  value.  This  effect  is  not  critical  in 
annual  energy  analysis  because  of  the  small  number  of  hours  for  which  it 
applies. 


1 

2 

3 

4 


TRNSYS  - A Transient  Simulation  Program  (Solar  Energy  Laboratory,  Uni- 
versity  of  Wisconsin,  March  1975). 

B.  Y.  H.  Liu  and  R.  C.  Jordan,  "The  Interrelationship  and  Character- 
istic Distribution  of  Direct,  Diffuse,  and  Total  Solar  Radiation," 
Solar  Energy,  Vol  4 (1960),  pp  1-19. 

A Test  Problem  and  Solutions  for  Solar  Heating  and  Cooling  Computer 
Simulation  Programs,  Draft  (U.S.  Department  of  Energy,  19/8). 

C.  M.  Randall  and  M.  E.  Whitson,  Jr.,  Hourly  Insolation  and  Mete- 
orological Data  Bases  Including  Improved  Direct  Insolation  Esti- 
mates, Aerospace  Report  ATR-78  (/592)-l  (The  Aerospace  Corporation, 
197?). 
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Figure  5.  Four  radiation  correlations. 


Time  of  Day  Factors 

Figures  7 and  8 show  August  radiation  data  separated  by  morning  and 
afternoon  observations.  The  H/HO  ratio  is  typically  higher  in  the  morn- 
ing, which  is  caused  by  a common  summer  weather  condition  in  Illinois. 
Mornings  start  clear,  but  small  cumulus  clouds  then  form  and  grow 
larger  in  size  and  number  throughout  the  day.  This  pattern  does  not 
occur  in  March  and  is  not  observed  in  the  data.  This  condition  can  be 
significant  for  summer  applications  of  beam  solar  radiation  and  will 
affect  the  orientation  of  the  solar  collector. 


Sources  of  Error 

The  severe  weather  in  February  affected  the  data  quality.  Long  pe- 
riods of  cloudy  weather  prevented  proper  alignment  of  the  shadow  band 
for  subsequent  clear  days,  especially  on  weekends  when  no  one  was 
present  to  adjust  the  band  as  soon  as  its  shadow  could  be  seen.  Occa- 
sionally, the  snow  made  it  too  hazardous  for  researchers  to  adjust  the 
shadow  band  on  the  roof. 

The  shadow  band  weathered  unexpectedly,  becoming  somewhat  reflec- 
tive instead  of  retaining  its  original  flat  black  color.  As  a result, 
it  now  reflects  radiation  from  the  ground,  especially  from  snow,  back 
onto  the  diffuse  measuring  pyranometer. 


MORNING  RRDIRTI0N : DRTS  2 


AFTERNOON  RADIATION:  OATS  215 


5 CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

1.  Insufficient  data  were  collected  to  establish  a definitive  cor- 
relation between  total  radiation  on  a horizontal  surface  and  its  direct 
and  diffuse  components. 

2.  The  close  agreement  between  the  correlation  developed  in  this 
study  and  the  SOLMET  correlation  indicates  that  the  former  should  be 
sufficiently  accurate  for  general  use  in  obtaining  information  about 
radiation  intensities  on  tilted  surfaces. 


Recommendations 

The  SOLMET  correlation  should  be  included  in  solar  simulation  pro- 
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